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FOREWORD 


This report documents source/radiation field coupling and workstation conversion for computer 
codes that are part of a fan tone noise prediction system based on linearized acoustic theory for 
turbofan and Advanced Ducted Propeller (ADP) engines. The study was conducted as Task X under 
the Aero-Propulsion Technology (APT) program, NASA Contract NAS3-25952. The NASA Project 
Manager for this contract is Mr. M. Murray Bailey, and the NASA Task Manager for Task X is 
Dennis Huff, NASA Lewis Research Center. Dr. D. C. Mathews served as Task Manager for Pratt 
& Whitney. 
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SECTION 1 


INTRODUCTION 

This document describes one portion of a fan noise prediction program designed to develop 
improved methodologies for source and noise field prediction. Other parts of the effort are 
described in separate reports. In particular, Reference 1 discusses improvements made to the 
Acoustic Radiation Code (herein referred to as the ARC), a finite element program for predicting 
noise radiation from a turbofan engine inlet. In Reference 2, evaluation of the ARC versus test data 
is described, along with work to update the wake model in the Bolt, Beranek and Newman/Pratt and 
Whitney Code (herein called the BBN/PWC). The BBN/PWC is a turbofan source noise code that 
was used to supply source input information to the ARC. Reference 3 presents work on a two- 
dimensional coupled cascade code that couples radiation between rotor and stator blade rows, while 
including swirl flow in the source region. Each of these efforts is aimed toward eventually obtaining 
better noise prediction tools for turbofan engines. The regions treated by the three codes are shown 
in Figure 1. 

This report concentrates on three areas: providing information so that developers of internal 
aerodynamics or aeroacoustic codes can understand the input requirements for the ARC, and also, 
the specific input formulation needed so that the BBN/PWC can be used as a source; development 
of a coupled interface between the source and radiation models; and modifications to the ARC for 
use on the Sun" and Silicon Graphics Iris" UNIX" workstations in use at Hamilton Standard, Pratt 
and Whitney, and NASA-Lewis. 

The reason for the first and last areas was to bring on board a versatile tool for studying noise. 
The ARC was a desirable choice because of its prior record of good correlation with test data and 
because its use of wave envelope elements reduces computational requirements. Since numerous 
source codes can be used to drive the ARC, general information is provided about the interface; 
however, since the BBN/PWC, in particular, is used as a source in work here and in Reference 2, 
more specific information is needed regarding its interfacing. The code conversion portion is 
discussed in the appendix. 

These efforts establish a noise prediction capability that progresses from the source region and 
then into the far field. Thus both the BBN/PWC and ARC regions shown in Figure 1 are now 
included in the predictions. The interface between the two regions is called the "fan face," in 
Reference 1, although the user may place this interface anywhere in the inlet. Hereafter, this cross 
section will be referred to as the "source input boundary.” 

The second area to be discussed is the coupling between the source and radiation regions. Work 
on it had its genesis in a desire to explore the possibility of obtaining more accurate predictions by 
feeding reflections back from the ARC to the source code. This work involved study of the waves 
reflected between the inlet back into the source region as well as development of a preliminary 
coupled code. This study was conducted in order to determine whether the contribution from the 
waves was sufficient to warrant the extra complexity of a coupled analysis. Related to this, there 
has been earlier work dealing with wave propagation and reflection in ducts (see refs. 4-8). This 
has been mainly analytical, providing reflection coefficients, but not indicating the effect on the 
source. The present work couples the ARC and BBN/PWC to obtain actual noise levels. 
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SUMMARY 


The Acoustic Radiation Code (ARC) is a finite element program used on the IBM" mainframe 
to predict far-field acoustic radiation from a turbofan engine inlet. In this report, requirements for 
developers of internal aerodynamic codes regarding use of their program output as input for the 
ARC are discussed. More specifically, the particular input needed from the Bolt, Beranek and 
Newman/Pratt and Whitney (turbofan source noise generation) Code (BBN/PWC) is described. In 
a separate analysis, a method of coupling the source and radiation models, that recognizes waves 
crossing the interface in both directions, has been derived. A preliminary version of the coupled 
code has been developed and used for initial evaluation of coupling issues. Results thus far have 
shown that reflection from the inlet is sufficient to indicate that full coupling of the source and 
radiation fields is needed for accurate noise predictions. Also, for this contract, the ARC has been 
modified for use on the Sun" and Silicon Graphics Iris" UNIX" workstations. Changes and additions 
involved in this effort are described in an appendix. 
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USING OUTPUT FROM SOURCE CODES 
IN THE ACOUSTIC RADIATION CODE 

In order for the ARC to predict the noise in the far field, it needs input information from a fan 
source code. To see what information might be useful to developers and users of source codes in 
providing such input, a number of internal unsteady aerodynamic methods (Euler, potential, 
acoustic) were surveyed. Some already had codes associated with them, while others could lead to 
engine source noise predictors in the future. These included work of Adamczyk, Mulac, and 
Celestina (ref. 14); Chi (ref. 15); Ni (ref. 16); Verdon, Barnett, Hall, and Ayer (ref. 17); and 
others. It was seen that output from source prediction methods must be processed or designed to 
be in the particular form required by the ARC. Therefore, general information is provided in the 
subsection below that should help developers of internal aerodynamic and acoustics codes understand 
the input requirements for the ARC. Additionally, because the BBN/PWC was used as a source 
for source/far-field coupling work at Hamilton Standard (see next section) and for code evaluation 
and noise prediction efforts at Pratt and Whitney (ref. 2), a specific procedure was needed for its 
input to the ARC. This information is presented in a second subsection. 


General Input Requirements 


The ARC requires source definition at an arbitrary plane in the inlet, the "source input 
boundary.” At this location, the pressure is considered to be the sum of a forward-going wave, P*, 
and a backward-going wave, P~. As mentioned previously, the user inputs mode information from 
P* and the ARC computes the radiation and the radial mode amplitudes for P~. The information 
provided here shows more explicitly what is required. Observe that, for later reference, both the 
forward- and backward-going waves are included in the discussion. However, only the forward- 
going wave is needed to provide input for the ARC. Note that not all source codes have pressure 
available explicitly and some may be two- rather than three-dimensional. It is the burden of the user 
to process the codes so that information can be extracted in the form required for the ARC. 


The full expression for the complex, nondimensional pressure, P, at the "source input boundary, " 
is given by 


P(r,4>,t) 


P*(r,4>,t) * P~(r,<t>,t) -SEE l p ~*m K« r ) CX P O' - W, - •tfO')) 

/•- a * *-1 

♦ P^* m («*/) exp - 7-u*. * *«lr)}] 


where the circumferential mode index, m, must satisfy the condition m = sB-kV (see the discussion 
regarding this in Section 4). In Equation 3.1, P is nondimensionalized using p m ci, where and 
c m are, respectively, reference values of the density and speed of sound. The cylindrical coordinates 
are shown in Figure 2. They are oriented as in the BBN/PWC. The quantity x 0 is the axial 
coordinate value at the "source input boundary," n is the radial mode index, s is the harmonic of 
blade passage frequency, B the number of rotor blades, V the number of stator vanes, N the total 
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SECTION 2 
BACKGROUND 

This section gives a very brief sketch of the two codes that have been used in this work, the 
ARC and the BBN/PWC. Complete details regarding the former can be found in References 1, 9, 
and 10; details for the latter are in References 11, 12, and 13. 

The ARC is a FORTRAN program used to study acoustic radiation from turbofan engines. It 
was developed at the University of Missouri-Rolla (UM-R) under the direction of Professor Walter 
Eversman. The code uses standard and wave-envelope type finite elements and an underlying 
potential mean flow to predict noise from axisymmetric ducts with or without centerbodies. 
Solutions calculated within the code are in terms of the velocity potential, $ m (x,r), associated with 
circumferential mode m, where the complete solution is $ = $„(x,r) exp { i (erf - Here x is 
the axial coordinate, r and are polar coordinates, u is rotor radian frequency, and t is time. 
Acoustic pressures are computed from <f>„. The program treats the ARC region indicated in Figure 
1. As input, both the index m and the reduced frequency (equivalently, rotor rpm or tip speed) must 
be set, along with complex forward propagating radial mode pressure amplitudes at the "source input 
boundary.” As mentioned in the previous section, the "source input boundary" is an arbitrary duct 
cross section where source information is input, and is the interface between source and radiation 
regions (see Figure 1). The ARC outputs complex backward propagating radial mode amplitudes 
at this same boundary. It also provides predictions for acoustic potential, pressure directivity, and 
acoustic pressure in the upstream radiation region. 

The BBN/PWC is a FORTRAN program that calculates the sound produced by the interaction 
of stator vanes with the mean velocity deficit wakes generated by the rotor blades of a turbofan 
engine. The wakes produce unsteady loading on the vanes, and this loading then serves as the 
source for the sound. The BBN/PWC was developed originally by Bolt, Beranek and Newman. 
Later, improvements were added by Pratt and Whitney. The fan/stator geometry for this study is 
shown in Figure 2. The code predicts the amplitudes of both upstream and downstream propagating 
modes, along with their sound power, in an infinite annular duct. To accomplish this, it uses the 
standard Green’s function approach. In this approach, the Green’s function is expressed as an 
expansion of normal duct modes; the loading (i.e., pressure distribution) exerted by the vanes on 
the fluid is evaluated for use in the Green’s function integral that gives the pressure anywhere in the 
duct. Knowledge of the upwash, i.e., the component of the wake velocity normal to the surface of 
the stator vanes, allows calculation of the pressure distribution on the vanes. This is done by 
dividing the stators into chord wise "strips" and computing the load on each "strip” assuming the 
stators to be af linear cascade of flat plates. By means of this process, once the upwash is known, 
the upstream and downstream pressure mode amplitudes are obtained in straightforward fashion. 
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where 


(3.7) 


r 2 * ('"*i («_/)«&• 


The coefficients P^, are not of concern to the user, as they are determined by the ARC as part of 
the solution. (Actually, it is the related velocity potential that is computed.) 


Using BBN/PWC Output 

When the BBN/PWC is used to generate thecomplex modal amplitudes, a straightforward 
conversion formula gives the modal information, to be used by the ARC. This is obtained 
by noting that, for a given circumferential mode order nr, the output pressure />£,, from the 
BBN/PWC, for a forward propagating wave at the "source input boundary,” is given by 

±r n - E A ~ ♦ b_ • < X8 > 

V Ip.lr >.1 


Observe that P£,, unlike in Equation 3.5, is dimensional; m and s denote the m* circumferential 
mode and the s* harmonic of blade passage frequency. In the above equation, p, is the nominal 
flow density and U is the mean axial fluid velocity. The square-bracketed terms represent 
expressions for the normal mode functions. They are the same as the in Equation 3. 1 except 
for the difference in the normalization for these between the BBN/PWC and the ARC. The 
coefficients A^, and are quantities that can be extracted from the BBN/PWC during a run. 
The values A^ are the upstream complex mode amplitudes, and the A^ and are the weighting 
constants for the normalized mode shapes. Note that for the BBN/PWC, the output that is used here 
is defined at the leading edge of the stator blade row. 

If one rewrites Equation 3.8 as 

*- * «_ • »•>_'))] < 3 - 9 > 

X npU a.| 


then the eigenmode in the square brackets is now normalized so that it has a maximum value of one. 
This is the salne as for the ARC. The quantity is the maximum value of the eigenmode over 
the range r„ ^ r ^ r®. 

For the ARC, the forward-propagating wave portion of Equation 3.5 can be rewritten as 

r~ - - E IT, - ^4 1’.f-w'-)!! (310) 

P-C. « 11(0 
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number of positive radial modes, 0 the rotor rotation rate, and t the time. The coefficients and 
f*imu 316 *h e complex pressure amplitudes for the modes and are nondimensionalized in the same 
manner as P. 

Further, the functions are the duct eigenmodes and are given by 

*»(*»/•) =r > - J J^L Y(tcr)} ( 3 . 2 ) 

*CO 


where J m and Y m are the standard Bessel functions of first and second kind of order m, is the 
eigenvalue for mode (m.ri), and J' m and Y' m are the derivatives of J m and Y m . The quantity T, is a 
normalization constant chosen so that the maximum value of ^fjr) over the range r H ^r<,r D is one. 
Here r 0 is the duct outer radius and r H is the hub radius. The parameters 7 ^ and 7 ^ are the 
forward-propagating and backward-propagating axial wave numbers for the (m.n)* mode and j* 
harmonic of blade passage frequency. They are given by 



r • 


- 

- 1 

MsBQ 


sB 0 

2 

- (i - m ! )«L 

l 

1 

M 

C ‘ \ 


* 4 


(3.3) 


where Af is the inlet axial flow Mach number and c a is the speed of sound in the duct. 

One can extract PJf), the complex amplitude of the (m.sf harmonic of P{r,^,t) with respect 
to r and f, by using the formula 

29 f 2 

w f W P(r,4>,t) exp {-/ (m<t> - sBQt)} d<t>dt (3.4) 

Once this is done, P M has the form 


P„(r) 


BQ 

(lx ) 1 


l 


N 


P„(r) - P^(r) . P„(r) - £ (3.5) 


*•1 


To drive the ARC with output from a source code not already written in modal form, the 
complex modal amplitude coefficients must be found from P^(r) by Fourier inversion. To 
accomplish this inversion, the user must first have separated out the forward-going portion, P* m , (r), 
of Pmj(r), using an approach appropriate for his code. Then the inversion formula is 

■ f M *„<«_') rdr (3.6) 
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SECTION 4 


DEVELOPMENT OF A COUPLED INTERFACE BETWEEN 
THE SOURCE AND RADIATION MODELS 

There are two approaches used in this report for linking source codes to the ARC. The first and 
simplest way is to input forward wave information directly from a source code to the ARC. This 
involves just one pass of each code and was the approach used at Pratt and Whitney for code 
evaluation (ref. 2). In this process, the ARC solution provides reflected wave information at the 
interface between the two codes, but the effect of the reflections on the source is ignored. This is 
the approach that was just described in the previous section. When the reflections are small, this 
approach should be satisfactory. 

Note that these reflected waves are the accumulated effect of reflection in the duct from the 
"source input boundary" to the inlet. From the mathematics, it is known that reflections will occur 
at any point in the duct where there are changes in geometry. These reflections have been studied 
in the past (see, for instance, refs. 4-8), but this report is believed to be the first to address the 
effects of the reflections on the source strength, i.e., blade loading. 

When reflections are large, a more sophisticated method might be preferable, one that would 
recognize waves passing the interface in both directions. This approach would iterate between 
source and radiation codes and, it is hoped, improve accuracy in both regions. 

As a first evaluation of reflection, the ARC alone was used. Reflections at the "source input 
boundary," due to unit modal inputs, were studied for a range of rotor tip speeds. Results from this 
investigation, described in the following subsection, demonstrated reflection large enough to justify 
proceeding further. 

A preliminary coupled code version was then developed by modifying the BBN/PWC to accept 
reflected waves and generate updated stator loading. The implementation of this approach is 
sketched out in the second subsection. Then the results of an initial set of exploratory runs, showing 
the effect of coupling on both the source and far-field regions, are presented in a last subsection. 


Reflection Ratio Curves 

Using onfy the ARC, curves of reflection ratios, | p"/p* | , versus engine tip speed, were 
derived (see Figures 3_through 6). Here, p and p+ are the reflected and incident complex pressure 
amplitudes, P^, and P^, at the "source input boundary." A unit modal amplitude was input for 
each choice of mode, (m,n), over a range of tip speeds. Reflected wave amplitudes for the first nine 
radial modes, for a given m, were obtained as output. 

Runs were for the fan inlet described by the finite element mesh in Figure 7, which is the figure 
used later to illustrate the quality of PostScript output. Dimensions shown are based on a unit fan 
tip radius. The model is for one of the test cases used by the UM-R in the past to test out ARC 
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if the normal mode functions, are written out as in Equation 3.2. 


By setting in Equation 3.9 equal to P^, in Equation 3.10, noting that the square-bracketed 
terms for each are equivalent, and equating coefficients of the bracketed terms because they are 
linearly independent functions, it follows that 


Jf 1 


- P.dPL 


(3.11) 


From this, the conversion relation becomes 



(f.m 



(3.12) 


or, using ARC and BBN/PWC program variables 

ZAI(n) * ( M ODIM)(PSIMAX)(AMNS) - 

(RHOFS)(CFSf ’ 

where ZAl(n) is RHOFS is p„, and CFS is c m in the ARC; MODIM is pJf/2, PSIMAX is 
and AMNS is in the BBN/PWC. The numerator of Equation 3. 13 is programmed into the 
BBN/PWC to output complex radial mode amplitudes in lb/in 2 (see ref. 2). This output can then 
be nondimensionalized using the denominator of Equation 3.13 to obtain ARC input This input 
information is supplied to the ARC via Input Card 4 of the ARC module PRATRADH. (See the 
appendix for a discussion of ARC modules.) 

Note that for coupling work at Hamilton Standard and for code evaluation work at Pratt and 
Whitney, the right-hand side of Equation 3. 13 was multiplied by 10,000 so that accuracy would not 
be lost in the ARC during calculations. Afterwards, outputs were scaled back by 0.0001 to produce 
actual values. 
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Table 1 

SCATTERING RATIOS FOR UNIT INPUT OF MODE (34,2) INTO (34,*) 



1,039 

1,069 

1,098 

1,106 

1,113 

1,128 

1,187 

1,247 | 

1 

.047 

.059 

E3 

.335 

.373 

.264 

.136 

1 

2 

.087 

.196 

■ 

1.991 

1.620 

.750 

.149 


3 

.011 

.012 

.042 

.046 

.050 

.038 

.117 

E3I 

4 

.001 

.005 

.001 

.003 

.004 

.004 

.011 

.020 

5 

.001 

.007 

.001 

.002 

.003 

.002 

.002 

.007 

6 

.001 

.001 


.003 

.003 

.002 

.002 

.003 

7 

.000 

.000 

.001 

.001 

.001 

.001 

.000 

.001 

8 

.001 

.001 

.001 

.002 

.002 

.001 

.001 

.002 

9 

.001 

.001 


.003 

.003 

.002 

.001 

.001 






*■2 cuton *>3 cuton 

(1,088 tUtec) (1,209 ft/aec) 


Table 2 

SCATTERING RATIOS FOR UNIT INPUT OF MODE (34,4) INTO (34,*) 
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development. All results were for an inlet Mach number, M, of 0.5, a free stream Mach number 
of 0.3, and a 38-blade ( B = 38) rotor. 

Figure 3 shows reflection ratios vs. V T for a fixed circumferential mode number of m = 25 and 
radial mode numbers of n - 1 to 4. Note that the peak of each curve occurs just a bit after the cut 
on tip speed of the corresponding mode. Reflection was considered to be important in regions 
where reflection ratios were 0.5 or greater. The horizontal bands in Figure 3 point out these 
regions. To quantify the widths of these bands, a bandwidth, W, was defined as 


W = 


AVJat \p-!p*\ =.5) 

— T - — — _ 1 lx 100% 

(; midband V T ) 


for each curve. In this expression, AV T (at | p'/p* | = .5) is the change in V r in passing from the 
left to the right side of the curves at a height of 0.5. As seen from Figure 3, these bandwidths vary 
from 21 % to 4% , where the bandwidth is greatest for n — 1. It can be seen that these bands cover 
more than half the speed range evaluated. We have looked at other m's and found similar behavior 
(see Figures 4 and 5). The curves in Figure 6 show the same point as above but for the situation 
where m varies and n = 1 is fixed. These curves, just as the previous ones, show that there is 
significant reflection in bands about cuton. Thus, Figures 3 through 6 show that there can be 
considerable reflection across the range of tip speeds important to turbofan engines. 


One might question points on the curves in Figures 3 through 6 where the ratios were greater 
than one. However, note that higher pressures in the reflected waves do not imply higher acoustic 
powers. Numerical experiments with a two-dimensional analysis indicate that, at M — 0.5, unit 
reflection of power produced roughly a factor of 2.5 for | p /p* | . 


Figures 3 through 6 deal only with reflection from a fixed radial mode n to n. There is also 
reflection from one radial mode to another. This is addressed in Tables 1 and 2 below. Data is 
from the same series of runs that produced Figures 3 to 6, and is for the case where m - 34. The 
entries in the tables are the ratios | p /p* | , where | p | now represents the reflected pressure 
coefficient for each of the nine radial modes indicated. The input, p*, for Table 1, is a unit 
amplitude for a relatively low order mode (34,2). For Table 2, the input is a unit amplitude for a 
slightly higher order mode (34,4). As can be sent from the data, there is scattering into both higher 
and lower order modes relative to the input mode. From the data in the tables, it is seen that there 
can be considerable reflection activity resulting from one radial mode reflecting into another. 

Observe that a stairstep type line has been added to both of the tables to indicate the division 
between modes that are cut on and cut off. Looking at ratios adjacent to these boundaries, note the 
tendency for the strongest reflected mode to move along these lines. This observation is based on 
only a small amount of data, so it is something we would wish to look into further in the future. 
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It can be shown that the upwash has the form (see ref. 11, Equation 82) 

w w = 53 w , ex P j**, z + - uflflr | (^-2) 

relative to stator vane v, v = 0, 1, . . . , V - 1, where s, B, V, and Q are the same as before. In 
Equation 4.2, an i has been added to the second term in the exponential; this was missing from 
Equation 82 of Reference 11. The w/s are the Fourier coefficients of the wake upwash at the j* 
harmonic of blade passage frequency, k, is the chord wise wave number, and z is the local coordinate 
along each of the stator vanes. The expressions for w s and k, are the same as in Reference 11 
except for differences produced by correcting for sign errors, by using the relative flow angle into 
the stator in place of the rotor stagger angle, and by shifting the z-origin from where the stator mid- 
chord meets the hub to where its leading edge intersects the hub. These modifications are needed 
to provide correspondence with the actual code. 

Using Equation 4.2, an expression can be obtained for Ap y , which is the s* harmonic of the 
pressure loading Ap appearing in Equation 4.1. It acts on vane v and (see ref. 11, Equation 86) is 
given by 

Ap, = P 0 U r w,f t (r, z) exp { -i 2t\ sB!V\ (4.3) 


where p„ is the nominal fluid density and U r is the fluid velocity relative to the stator vanes. The 
minus sign in the exponential corrects a sign error in Reference 11. The quantity f,(r,z) is the 
elemental stator vane loading function which comes from solving 

exp{i*,z} = ^ K e (z-y)f,(r,z)Q (4.4) 


Here K t {z-y) is the kernel function for a linear cascade and is discussed in Reference 11, 
Appendix B. The integration is from 0 to 2b, rather than from -b to b as in Reference 11, because 
of the shift in origin for z mentioned previously. 


By next returning to Equation 4.1, expanding G there in terms of eigenmodes, and then using 
Fourier transform relations, one can obtain an expression for the complex pressure modal 
amplitudes, pZ*,. This expression will depend on A p v above. Using Equation 4.3 in this expression, 
one then obtains the relation (see ref. 11, Equation 187) 


<r c (U,fV) 2 V ,x 


f 1 

W, 

J 

U r 


r 


*»(*«.*)(— «> s «, + sino 0 

o 


exp | -i -^5 2 j <t c j BC^i^da 


(4.5) 


for the amplitudes, pZu, nondimensionalized by 'hpJJ 1 . The subscripts m, n, and s signify the 
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As a final note, observe that the shape of the curves in Figures 3 through 6 can be explained as 
follows. An input wave can be shown to propagate axially with complex amplitude proportional to 

the term exp{-w,Kjjc - . Here A, and A 2 are positive constants, V TCO is cut on tip 

speed, and x is the duct axial coordinate. This expression comes from rewriting the factor 
exp {-/ in Equation 3.1; x a w as generali zed here to x and was rewritten in terms of V T 

and Vt.co- Just below cuton, the *} portion of this exponential is real and nearly one. 

Hence, one would expect the input wave to propagate through the duct, from the "source input 
boundary" to the inlet, nearly undiminished, and one would expect the reflected wave to be largest 
in this tip speed region. Because this exponential decreased rapidly for V T just a bit below cuton, 
the drop-off to the left of V TCO should be steep, which it is. 

To the right of V TCO , the square root part of the exponential becomes imaginary. Thus the whole 
exponential is sinusoidal and the input wave propagates toward the inlet undiminished. However, 
as V T increases, the axial wavelength that is proportional to 1 / (*,K r - >^») decreases relative 

to die annular width of the duct. As the wavelength decreases, more of the input wave will be 
radiated outside the inlet and less will be reflected. Hence the ratio, | p'/p* | will drop off to the 
right of V TCO , as is seen in the curves. 


Implementation of Coupling 

To implement coupling, the source code was modified to make use of the reflected wave 
information from the ARC to update loading on the stator vanes and, from this, the modal 
amplitudes. The source code chosen for this purpose was the BBN/PWC. Because it was being 
used at Pratt and Whitney, it could be quickly brought on board; also, because it was based on the 
use of duct modes, as is the ARC, the mathematics for coupling was highly compatible. Before 
briefly describing how this was done, the approach used for the original derivation will be 
summarized here. The reader is referred to Reference 11 for further background. Note that the 
coupling implemented here is used for the runs described in the next subsection. These runs 

comprise the first step of an iteration process and are explored to establish the importance of 
coupling. 


As mentioned in Section 2, the original BBN/PWC predicts the amplitudes, of the pressure 
modes that are generated by the loading on the stator blades coming from the upwash, vF, of the 
rotor wake. One starts with the Green’s function representation for pressure (see ref. 11, Equation 
4.1), which is given by 


P(x,t) » - J * | *00 • VG(x,y,r - r) Ap (y,r) dS(y)dr 


(4.1) 


where S is the stator vane surface area, x and y are coordinate vectors (in cylindrical coordinates) 
in the duct, n is the unit normal to the stator blades, G is the Great’s function, V is the gradient 
operator, / and r are time, and &p is the pressure loading on the vanes. 
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The velocity, v*, calculated from is given by v* = V (r,<t>jc). This gives an expression 
for the velocity anywhere in the duct. If the component of v* normal to the stator vanes is evaluated 
along the stator vanes, one has w*. After proceeding through this calculation, coordinates and 
variables must be converted to those used in the BBN/PWC so that results can be used in this code. 
In doing this, differences in coordinate orientation, scaling, and definition of similar variables must 
be taken into consideration. Also, correct phasing must be maintained. The result is the relation 

w R = ex P z ~ i “ isBQt | 

for the upwash, due to reflected waves, on vane v, v = 0, 1, . . . , V - 1. 

In Equation 4.9, the w/s are the Fourier coefficients of the reflection upwash at the harmonic 
of the blade passage frequency for radial mode./, and kp is the chordwise wave number. The 
quantities and k p correspond to w, and k, in Equation 4.2, but the expressions defining these are 
obviously very different. The quantities above depend on r, a„ 7^», and the vane sweep 
parameters. 


If one proceeds from this point as previously for the wake case, one then obtains the pressure 
modal amplitudes, /£,, representing the effect of reflections from the inlet striking the stator and 
re-radiating upstream. They are nondimensionalized by l hp e lP and are as given below: 



'Ap.U 2 


Q'iVJUfV j, r* 

t (1 J- 1 l. 





* 




(4.10) 


Just as for the wake case, the mathematics that gives Equation 4.10 requires that m be such that 
m - sB - kV, where k may be any integer. The functions C^/o) depend on, among other things, 
the geometry, wave numbers, and the elemental vane chordwise pressure distributions, ( r,z ). 

Analogous to the situation for the/ (r,z) in the wake case, the/, (r,z) are calculated by solving the 
integral equation 

«p{»v> * ^ (411) 

where K (z-y) is the kernel function employed before. Observe that there are no phase terms in 
Equation 4. 10 similar to the one in Equation 4.5. This is because such terms have been absorbed 

into the w^’s in Equation 4. 10. 

As a last comment, note that the modifications described above have been incorporated into the 
SRCHV, FILLOV, PRES, OPFILE, and CLFILE subroutines of the BBN/PWC. 
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(m,n) mode at the j* harmonic of blade passage frequency. In Equation 4.5, U is the axial fluid 
velocity; a c = 2 bj/r D where b T is the vane semi-chord at the tip and r D is the duct outer radius; a r 
= r n/ r o is toe hub radius; and 7 «« = y^o where y^, the axial wave number, is the same as 7 * r 
in Equation 3.3. In defining y^, 7 ^ or y^, is picked according to whether is for an 
upstream-going wave or for a downstream-going wave. Further, in Equation 4 . 5 , = * r 

where K", is the square root portion of 7)W ; is the normal mode; where is the 

eigenvalue for mode (m,n); a = r/r^ a, is the stator stagger angle; 5 t and 6 2 are 

nondimensionalized vane sweeps; B = b/b^ where b is the vane semi-chord; and the are the 

chordwise functions appearing in Reference 11, but with minor modifications. Note that pi, in 
Equation 4.5 is nondimensionalized using Vip.l/ 2 , rather than p g lf as in Reference 11, so thaTthe 
expression corresponds to that in the code; also, typographical errors have been corrected. Lastly, 
observe that the mathematics to arrive at Equation 4.5 imposes the condition that m = sB - kV, 

where k may be any integer. The right-hand side is opposite in sign to the condition in Reference 
11 , because of an error there. 

For the modified version of the BBN/PWC, an additional upwash, 1 /, due to reflected waves, 
must be calculated. The noise is now produced by an upwash, w, on the stator, given by 

w = w w + w* ( 4 - 6 ) 

The upwash w w is the same as before; the additional upwash, w*, is calculated using the reflected 
wave information provided by the ARC. Because each ARC run is for a fixed frequency, u = sBQ, 
w* will not be given by a sum of terms over s as was (see Equation 4.2), but will be given by 
only one such term. Once w* is found, the pressure amplitude coefficients, pi,, can be determined 
by a process parallel to that used above to derive the pi,. Then the total pressure amplitude is 

Pmtu ~ Puna + Pmtu ( 4 . 7 ) 

To obtain w* , one starts with the velocity potential for the reflected wave, which is given by 

N 

(kjr ) exp {/ {-nub -yL* + sBQr)} (4.8) 

>-■ 

The variables and coordinates here are those used by the ARC. This means that the orientation for 
( r ’4>’ x ) and Q is opposite to that shown in Figure 2, which gives the orientation used by the 
BBN/PWC. Variables have the scaling used by the ARC; thex-origin, the axial position where the 
"source boundary interface" is located, falls where the leading edge of the stator intersects the hub. 

The subscripts m and s, on the left, are the indices assumed for the ARC run. The coefficients 
are the velocity potential complex modal amplitudes from the solution for the reflected waves. 
There are N of these, corresponding to the number of positive radial modes, indexed by j, chosen 
for the ARC run. The modes can be cut on or cut off because both types are permitted by the 
ARC. However, the BBN/PWC only computes cut on modes. Hence, the loading generated by 
Equation 4.8 will only contribute to those radial modes, specified below by the index n, which are 
propagating modes in the BBN/PWC. 
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somewhat removed from cuton, the patterns are essentially the same but the levels are different. 
At the second tip speed, nearer to cuton, both the shapes and the amplitudes vary. 

Current efforts show that coupling is important. However, present results are for only one 
iteration step. It will be important in future work to continue the iteration and study convergence. 
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Coupled Runs 


The calculations presented in the first subsection showed there is considerable reflection activity 
in the duct. To further explore the importance of this, an initial study of the effect of coupling on 
noise levels was begun. To this end, a series of computer runs was made that coupled the 
BBN/PWC with the ARC. For a given tip speed (actually rpm in the BBN/PWC and reduced 
frequency in the ARC), each run involved first executing the BBN/PWC, then feeding its output as 
input to an initial ARC run. Reflected wave output from this ARC run, at the "source input 
boundary, was next fed into a second BBN/PWC run where it provided updated stator blade 
loading. The resulting BBN/PWC output was then used as input for a second ARC run. 

These initial tests were made for the same geometry (see Figure 8) as used in the evaluation 
studies at Pratt and Whitney (ref. 2) for a similar tip speed range, 9,300-12,000 ft/sec. The setup 
treated a 16-blade, 22-vane, two-chord spacing, medium length ADP inlet model. The exterior 
Mach number was 0.2. The circumferential mode number studied was m = -6, also studied by 
Pratt and Whitney, and only the first BPF harmonic was used. Essentially the same input as used 
at Pratt and Whitney for the V T = 788 ft/sec case was used as a starting point. The rest of the 
coupled runs employed the same properties, changing only the tip speeds. In the coupled 
BBN/PWC runs, the first nine radial modes for tn — -6, not all of which were cut on, were used. 
This is permissible (see the discussion in the previous subsection) even though the BBN/PWC only 
gives output for cut on modes. 

Results are shown in Figures 9 through 12. Figure 9 compares total upstream sound power in 
the first and second BBN/PWC runs, i.e. , before and after coupling, as a function of tip speed. 
Note that cuton for the second radial mode, (-6,2), is at 706 ft/sec where the two curves peak. 
Below this speed only radial mode (-6,1) contributes to power; above this point both modes (-6,1) 
and (-6,2) contribute. Observe that radial modes (-6,1) and (-6,2) here correspond to (-6,0) and (- 
6,1) in Reference 2. It is seen that in a band of tip speeds above second radial mode cuton, there 
is a sizable increase in total power from the uncoupled to the coupled case. (Note that the increase 
just below cuton of the (-6,2) mode here and in Figure 10 should be disregarded; we trust that had 
there been additional points computed in this vicinity, there would have been vir tuall y no delta.) 
At a typical point, V T = 788 ft/sec, a bit removed from cuton, the increase is 5.7 dB. 

Figure 10 shows a similar comparison for the ARC between the first and second runs, where the 
first is before coupling and the second is after. The ordinate here is maximum sound pressure level 
at 10 inlet radii, rather than total power as before. This quantity is used because it, rather than total 
power, is easily derivable from ARC runs. Also, it was felt that this parameter gave a good 
indication as tb whether coupling was, or was not, having an effect Just as for the prior BBN/PWC 
case, plots here show a band above cuton where the change in noise levels due to coupling is 
significant, and at V T = 788 ft/sec, the increase in SPL was comparable to the power increase 
computed via the BBN/PWC. 

Figures 11 and 12 illustrate further the effect of coupling on far-field noise. Figure 11 shows 
the change in pressure directivity patterns in passing from the first to the second ARC run at 
V T = 788 ft/sec. Figure 12 gives the same information at V T = 710 ft/sec, which is nearer to 
cuton. Both show that coupling has made a major difference in noise results. At the first tip speed, 
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SECTION 6 


APPENDIX 

CONVERSION OF THE ACOUSTIC 
RADIATION CODE TO THE WORKSTATION 

This appendix discusses modifications that were made to the ARC so that it would run on Sun 
and Silicon Graphics Iris UNIX workstations. The changes were of two types. The first set 
consisted of those that were needed to bridge the gap between what was acceptable for the 
UNIX/FORTRAN workstation-based environment as opposed to the IBM mainframe one, which was 
more forgiving. The second, and main change, involved adding subroutines to provide the 
PostScript" plot output requested by NASA. The original code provided CalComp’ output. 

The workstation version consists of the five modules listed below. The designations for them 
are the same as for the UM-R modules except that an "H," for Hamilton Standard, has been added 
to the nam es of the workstation versions to distinguish between the two. 

PRATREH - generates input data for nacelle and 

centerbody geometry for use in PRATMESHH 

PRATMESHH - generates the finite element mesh and 
the duct eigenvalue solution 

PRATFLOWH - obtains the potential flow solutions 

PRATVELH - superposes the solutions above for 
use in PRATRADH 

PRATRADH - generates the acoustic radiation solutions 

Note that there are two UM-R modules, PRATRADA and PRATRADB, corresponding to the 
workstation version, PRATRADH. The first of these uses CalComp plot routines, while the second 
creates plot files for a commercial program called TECPLOT. Because a non-proprietary approach 
was preferred, PRATRADA, rather than PRATRADB, was chosen for workstation conversion. 
Also note that the input and output descriptions for both the workstation and IBM mainframe 
versions of thf above modules are the same. See Reference 1, Appendix A, for a description of the 
input and output for the modules. 


This subsection discus*** changes that were needed because some coding, though acceptable in 
the IBM main frame environment, would not work in the Sun/Silicon Graphics one. Of the five 
modules above, PRATPREH remained essentially the same as the original. The other four required 
changes. 
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SECTION 5 

CONCLUDING REMARKS 

In this report, information was presented regarding use of output from source codes as input for 
the ARC (Acoustic Radiation Code). This involved describing general input requirements for any 
source code as well as for specific ones for the Bolt, Beranek and Newman/Pratt and Whitney Code 
(BBN/PWC). A result of this effort is a one-pass capability using both the BBN/PWC and the ARC 
to predict noise in both the source and far-field regions. Development of a coupled interface 
between the source and radiation models was described. In this regard, the BBN/PWC was 
modified to accept reflected wave input from the ARC. This input was used to update blade 
loading, by means of which noise predictions were updated. In the appendix which follows this 
section, conversion of the ARC from the University of Missouri-Rolla IBM mainframe version to 
the Sun and Silicon Graphics Iris UNIX workstation versions is discussed. 

The coupling process and reflection in the duct were studied. Using the ARC alone, unit 
amplitude forward propagating modes were input at the "source input boundary." The code 
computed the reflected amplitudes for the same radial modes as those input. The ratio of these 
amplitudes, | p‘/p + ] , versus rotor tip speed, for a selection of circumferential mode and radial 
mode pairs, (m,n), indicates that in significantly wide bands about cut on tip speed, there is 
significant reflection from the inlet. There was also significant reflection into other radial modes 
in the vicinity of cut on boundaries 

Additional runs were made using the BBN/PWC-ARC preliminary coupled version of the code 
for an ADP model, with circumferential mode order m — - 6. For this process, a first BBN/PWC 
run was made. Output from this first run was fed into the ARC. Then this BBN/PWC-ARC cycle 
was repeated. Plots of total sound power versus tip speed for propagating waves in the duct and 
of maximum sound pressure level versus tip speed for waves at 10 inlet radii in the far field, 
showed significant increases in radiation levels in a band above second radial mode cuton. For 
instance, at a tip speed of 788 ft/sec (cuton was at 706 ft/sec), the increase in total power in the duct 
was 5.7 dB; a similar increase in SPL was also seen in the far field. 

Based on results thus far, we can conclude that coupling reflection back to the source does have 
an effect on net radiation and should be included in predictions. Note, however, that work here has 
been for only one iteration step. In the future, further iterations need to be studied and convergence 
should be investigated. 

In view of* the above remarks, there are a number of recommendations re gardin g further work 
in this area. Coupling should continue to be studied by further developing the code that links the 
ARC and the BBN/PWC. Changes in stator loading, sound pressure, or power in the source region 
and far field, and scattering between radial modes should be investigated further. By means of 
additional runs over a range of tip speeds of interest to ADP engines, additional light should be shed 
on physical, computational, and mathematical aspects of coupling. As part of this additional work, 
behavior close to additional cut off points should be studied and, as pointed out above, convergence 
issues should be considered. 
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The second difference does not involve code. The Silicon Graphics compiler interprets record 
lengths of unformatted direct access files as numbers of words rather than numbers of bytes, as for 
the Sun. Code could have been changed for the Silicon Graphics to account for this; however, to 
keep both versions as nearly the same as possible, this was not done. By using the "-old_rl" option 
during compilation, record lengths are interpreted as bytes rather than words, thus accomplishing 
the same end. 


Adding PostScript Output Capability 

Plotting in the original code was accomplished through calls to standard CalComp routines. To 
provide PostScript plots, instead, PostScript subroutines were written to take the place of the 
CalComp equivalents. These generate a plotting file, fort. 14, for each of the modules, 
PRATMESHH, PRATFLOWH, PRATVELH, PRATRADH, that later can be plotted out on a 
PostScript plotter or viewed on a monitor by means of one of the utility programs available for the 
Sun or Silicon Graphics workstations. Examples of PostScript plots are shown in Figure 7. 

The PostScript subroutines that were added are listed below. They have the same arguments and 
perform the same functions as the original CalComp subroutines. 

SYMBOL (X, Y, HEIGHT, IBCD, ANGLE, NCHAR) 

- prints text or special symbols 

PLOT (X, Y, IPEN) 

- moves pen, either up or down, in a straight 
line to a new location 

PLOTS (I, J,K) 

- initializes the PLOT subroutine 

NUMBER (X, Y, HEIGHT, FPN, ANGLE, NDEC) 

- prints a number 

Arguments for the above subroutines are defined as follows: 

X, Y - coor dinates to which pen is moved or 
where symbol or number is placed 

HEIGHT - height of character or number 

IBCD - the text or integer to be printed 

ANGLE - angle from the horizontal at which to print 
text or number 

NCHAR - number of characters to be printed 
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The most changes were needed to correct difficulties with plotting. In this connection, some 
general changes related to maintaining consistent data-typing were required. These were to insure 
that particular variables were of the same data-type and precision, double or single, when used in 
different routines, and also to make sure that there were no conflicts between the arguments in the 
calling statements for subroutines and the argument definitions in the subroutines themselves. 

In particular, "IMPLICIT REAL-8 (A-H, O-Z )" and "IMPLICIT INTEGERS (I-N)" 
declarations were added at the beginning of most routines and function subprograms. Also all 
REAL, INTEGER, and COMPLEX" data-type declarations for particular variables were 
standardized to be "REAL*8," "INTEGER-4," and "COMPLEX-16." Additionally, double- 
precision constants were written in exponential "D-notation" (e.g., 12.5D+0) in subroutine 
argument lists, function calls, and elsewhere. 

Further adjustments were also needed to eliminate problems in plotting. In PRATFLOWH, 
PRATVELH, and PRATRADH, the subroutine CONTUR needed changes. The same was true for 
the subroutine CONTP in PRATRADH. Both of these subroutines assist in contour plotting and 
were causing nondesired coordinates to be defined and plotted together with intended ones when they 
called the subroutine NWPLT2. Note that NWPLT2 is used to plot different types of curves (solid, 
chain-dotted, and dotted). Difficulties were cleared up by changing "C ALL NWPLT2 (XP RP 
2, ITYPE, PLOTD)" to "CALL NWPLT2 (XP, RP, 3, ITYPE, PLOTD)" in CONTUR and 
CONTP, and changing the coordinates XP(2) and RP(2) to XP(3) and RP(3). "INPT=2" was 
changed to "INPT=3." In NWPLT2, the statement "NS=NP-IDBL" was modified to be 
"NS = NP"; NS was a limit index-value for a DO-LOOP where coordinates were defined. This 
latter change was also made in the subroutine NWSCLE which was used in all the ARC modules 
except PRATPREH; NWSCLE sets up suitable ranges for a given set of points when plotting. The 
purpose of this change was to be sure that coordinates were properly defined in NWSCLE also. 

Regarding nonplotting-related changes, in the EIGEN subroutine of PRATMESHH, there was 
originally a call to the IBM mainframe service routine, XU FLOW. This suppresses interrupts to 
the program from exponent underflow, which could be disruptive. The call to XUFLOW was 
disabled in the workstation version because it is not available there. Nothing was added to replace 
this because default 1 EEF . floating point exception handling for UNIX workstations does nearly the 
same thing as XUFLOW. 

There were a number of other, generally minor, changes made to the code, even though they 
would have no immediate effect on workstation runs. The first of these was implemented with the 
thought of eliminating possible future difficulties should the code be moved, later, to other computer 
platforms. Most of these dealt with eliminating additional minor data-typing inconsistencies between 
calling programs and subroutines. These were suggested through tests using the Lahey" FORTRAN 
compiler on the IBM PC" and using the ICA (Intercompilation Analysis) option during IBM 
mainframe compilations. 

Note that there are two slight differences between the Sun and Silicon Graphics versions of the 
program. The first regards coding where, for the Sun, in PRATRADH, an "OPEN (UNIT =6, 
FORM = ‘PRINT’)" statement was added to prevent unwanted control characters from being printed 
out in hardcopy. The option "PRINT" is not available for Silicon Graphics workstations. 
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I 


IPEN 

- 2 = pen down; 3 = pen up; -2 
establishes new origin 

I,J 

- no longer used; set to 0 

K 

- logical output device number 

FPN 

- floating point number to be printed 

NDEC 

- controls precision of number printed 


Further details regarding the functions and arguments of the above subroutines can be found in the 
CalComp manual (ref. 18). 

Note that in transferring to PostScript from CalComp, it was necessary, also, to modify a very 
small amount of plotting-related code in some of the ARC modules. This was to adjust the spacing 
for text and numbers so that they would fall in the proper locations under PostScript. Keeping the 
original CalComp placement would not have done this. 
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LIST OF SYMBOLS 


b 

b 

b T 

Co 

Co. 

exp{ } 


fj. (r,z) 

i 

K 

K 



m 

n 

n 

P 

P' 

P* 


vane semi-chord 
b/b T 

vane semi-chord at tip 

speed of sound in duct 

ARC reference value for speed of sound 

exponential function 

elemental vane chordwise pressure distribution for s® harmonic of blade passage 
frequency 

elemental vane chordwise pressure distribution for/ 1 radial mode and s® harmonic 
of blade passage frequency 

vane chordwise wave number for s® harmonic of blade passage frequency 

vane chordwise wave number for y® radial mode and j® harmonic of blade passage 
frequency 

square root portion of 
K~r D 

circumferential mode index 
radial mode index 

unit vector in direction normal to stator vanes 
acoustic pressure 

complex amplitude for reflected pressure wave at "source input boundary" 
complex amplitude for incident pressure wave at "source input boundary" 
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A., 

AMNS 

B 

BMN 

CFS 
C„{ ) 
Cj^i ) 
G( ) 

J m { ) 

j: 

k c ( ) 

M 

MODIM 

N 

P(r.4>,Q 

P~(r,<t>,t) 

P * (r,4>,t) 


first weighting constant used in BBN/PWC to normalize eigenmode (m,n) 

BBN/PWC program variable for nondimensional upstream pressure complex 
amplitude for mode (m,n) and s* harmonic of blade passage frequency 

nondimensional upstream pressure complex amplitude for mode (m,ri) and s* 
harmonic of blade passage frequency in BBN/PWC 

number of rotor blades 

BBN/PWC program variable for second weighting constant used in normalizing 
eigenmode (m,n) 

second weighting constant used in BBN/PWC to normalize eigenmode (m,n) 

ARC program variable for c„ 
function used in Equation 4.5 
function used in Equation 4.10 

time-dependent Green’s function for hard-walled annular duct 

Bessel function of the first kind of order m 

first derivative of J m 

kernel function for linear cascade 

inlet axial flow Mach number 

BBN/PWC program variable for p,l/*/2 

total number of positive radial modes for ARC run 

ARC complex nondimensional pressure at "source input boundary" 

pomplex nondimensional pressure for backward-propagating wave at "source input 
boundary" 

complex nondimensional pressure for forward-propagating wave at "source input 
boundary" 
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complex amplitude of acoustic pressure, produced by reflected waves, for mode 
(m,n) and 5 ttl harmonic of blade passage frequency 

complex amplitude of acoustic pressure, produced by wake, for mode (m,n) and 
harmonic of blade passage frequency 

polar radius in cylindrical coordinates 

duct outer radius 

hub radius 

index for harmonic of blade passage frequency 
time 

stator vane index number 

velocity in duct produced by reflected waves 
w w + w * 

Fourier coefficient of wake upwash on stator vane for s* harmonic of blade 
passage frequency 

upwash along stator vanes produced by reflected waves 
upwash along stator vanes produced by wake 

Fourier coefficient of reflected wave upwash on stator vane for j* radial mode and 
r* harmonic of blade passage frequency 

axial coordinate in cylindrical coordinates 

field point for Green's function in cylindrical coordinates 

axial coordinate value at the "source input boundary” 

vane chordwise variable 

integration point for Green’s function in cylindrical coordinates 
vane chordwise variable 



Y m ( ) 

£( ) 
w 

ZAl(n ) 




7, 


'Yjmu 


Tmw 


Ynuii 


r, 

r 2 


5i, 5j 


V 

Pa 

Pa. 

a 

°c 


Bessel function of the second kind of order m 
first derivative of Y m ( ) 
bandwidth for reflection ratio curves 

complex amplitude for forward nondimensional pressure wave, with radial mode 
number n, at the "source input boundary" 

vane stagger angle 

axial wave number for (/n,n)* mode and 5* harmonic of blade passage frequency 

1*uJd 

axial wave number for reflected wave for mode (m,n) and s* harmonic of blade 
passage frequency 

axial wave number of forward-propagating wave for mode (m,ri) and s* harmonic 
of blade passage frequency 

normalization constant used in Equation 3.2 

normalization constant defined in Equation 3.7 

pressure loading on stator vanes in time domain 

pressure loading on stator vane v in frequency domain 

nondimen sionalized stator vane sweeps 

gradient operator 

eigenvalue for mode (m,n) 

BBN/PWC nominal flow density 

ARC reference value for density 

r/r D 

2bj/r D 


- 28 - 



P„(r) 


P~ m (r) 

PL(r) 


P- 


] 


complex nondimensional amplitude for (m.s)* harmonic of P(r,<t>,t) with respect to 
r and t at "source input boundary 


backward-propagating wave portion of P m (r) at "source input boundary" 

forward-propagating wave portion of P M (r) at "source input boundary" 

complex amplitude for harmonic of P (r,<f>,t) with respect to r and t for 
forward-propagating wave at "source input boundary” 

reflected wave complex nondimensional pressure amplitude for mode (m,n) and j* 
harmonic of blade passage frequency, as used in F/ p»a H on 3.1 


complex nondimensional pressure amplitude for forward-propagating pressure wave 

for mode (m.n) and j* harmonic of blade passage frequency, as used in Equation 
3. 1 


complex nondimensional pressure amplitude for reflected wave for mode (m.n) and 
s* harmonic of blade passage frequency, as used in Equation 3.5 

complex nondimensional pressure amplitude for forward-propagating pressure wave 
for mode (m.n) and s* harmonic of blade passage frequency, as used in Equation 


PSIMAX BBN/PWC program variable for 

RHOFS ARC program variable for p„ 


S 


stator vane surface area 


U BBN/PWC mean axial fluid velocity 

U r fluid velocity relative to stator vanes 

V number of stator vanes 


V T rotor tip speed 

Vt.co rotor cut on tip speed 


- 27 - 



ACKNOWLEDGEMENT 


The author would like to thank D. B. Hanson for suggesting the project treated here and for 
his support and helpful suggestions. 


- 30 - 



o T rJr D 

t time 

4> polar angle in cylindrical coordinates 

* ARC velocity potential, with space and time dependency 

(x,r) (x,r)-dependent portion of $ 

ARC velocity potential, with space and time dependency, for reflected waves 

complex amplitude for ARC reflected wave velocity potential for mode (m,n) and 
s* harmonic of blade passage frequency 

*„( ) ARC normal duct mode for circumferential mode m 

¥*« maximum value of the duct eigenmode over the range r„ £ r £ r D for BBN/PWC 

u radian frequency 
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Figure 1. Fan noise design/prediction sy 
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(ref. 2). 
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Figure 4. Reflection ratios. 
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Figure 3. Reflection ratios. 
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Figure 6. Reflection ratios, m varies 



Figure 5. Reflection ratios, m = 34, n = 1, 2, and 4. 






Medium inlet 



Figure 8. Advanced ducted propeller inlet configuration (ref. 19). 
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Figure 7. Typical PostScript plots. Finite element mesh within the inlet (nomenclature 
added) and polar SPL directivity. 
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Figure 10. Comparison for ARC runs. Sound Pressure Level at 10 inlet radii for m - 
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Figure 9. Comparison of BBN/PWC runs. Total upstream power for m = 




(19.74 x 10^ 
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